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Background: Atomic nuclei often exhibit collective rotational-like behavior in highly excited states, well above 
the particle emission threshold. What determines the existence of collective motion in the continuum region, is 
not fully understood. 

Purpose: In this work, by studying the collective rotation of the positive-parity deformed configurations of the 
one-neutron halo nucleus ^^Be, we assess different mechanisms that stabilize collective behavior beyond the limits 
of particle stability. 

Method: To solve a particle-plus-core problem, we employ a non-adiabatic coupled-channel formalism and the 
Berggren single-particle ensemble, which explicitly contains bound states, narrow resonances, and the scattering 
continuum. We study the valence-neutron density in the intrinsic rotor frame to assess the validity of the adiabatic 
approach as the excitation energy increases. 

Results: We demonstrate that collective rotation of the ground band of ^^Be is stabilized by (i) the fact that 
the £ = 0 one-neutron decay channel is closed, and (ii) the angular momentum alignment, which increases the 
parentage of high-£ components at high spins; both effects act in concert to decrease decay widths of ground-state 
band members. This is not the case for higher-lying states of ^^Be, where the i = 0 neutron-decay channel is open 
and often dominates. 

Conclusion: We demonstrate that long-lived collective states can exist at high excitation energy in weakly bound 
neutron drip-line nuclei such as ^^Be. 

PACS numbers: 21.60.Cs, 21.10.Re, 21.10.Tg, 24.10.Eq 


Introduction. - Studies of exotic nuclei far from the 
valley of beta-stability reveal novel features, such as the 
formation of halo structures [1, 2], near-threshold cluster¬ 
ing effects [3-6], and presence of new types of correlations 
[7]. In all these cases, the atomic nucleus exhibits prop¬ 
erties characteristic of an open quantum system, whose 
properties are dramatically altered by the coupling to 
scattering and reaction channels [8]. 

While the impact of reaction channels on nuclear struc¬ 
ture has been recognized [9], the fact that some highly 
excited nuclear excitations can be interpreted in terms of 
nuclear clusters and nuclear molecules [3, 4] that experi¬ 
ence collective motions such as rotations and vibrations, 
is quite astonishing. The success of the collective model 
of atomic molecules relies on the validity of the adiabatic 
Born-Oppenheimer approximation [10] that relies on the 
vast difference between time scales for single-electron and 
ionic motions. However, when it comes to atomic nuclei, 
the time separation between single-nucleon and collec¬ 
tive motion is small, and the adiabatic approximation is 
expected to be badly violated [11]. How come, there¬ 
fore, that a highly excited state, undergoing rapid par¬ 
ticle emission, can be viewed in simple geometric or al¬ 
gebraic terms involving rotating or vibrating fields, or 
potentials, common for all nucleons? An excellent re¬ 
cent example is offered by the spectrum of discussed 


in the geometric/algebraic language of three a-particle 
arrangements [12]. 

Whether a broad feature observed in scattering exper¬ 
iments can be understood in terms of a nucleus, or a nu¬ 
clear state, is a subject to ongoing debate. Consider the 
single-particle time scale. The average time it takes a nu¬ 
cleon to go across a light nucleus {A « 10) and come back 
can be roughly estimated at Tg p, « 1.3-10"^^ sec [13], and 
corresponds to the time scale needed to create the nuclear 
mean field. If one relates the half-life of the unbound nu¬ 
clear state to its decay width T through T 1/2 = ln(2)?i/r, 
one is tempted to conclude that broad scattering features 
with Tif 2 < Tg.p. (or T > 3.5 MeV for A « 10) can hardly 
be interpreted in terms of nuclear states [14]. Of course, 
there is no sharp borderline that separates bona fide nu¬ 
clear states from broad features seen in scattering exper¬ 
iments, and this often results in interpretational difficul¬ 
ties [15]. To cloud the issue even more, a quantitative 
experimental characterization of broad resonances, em¬ 
bedded in a large non-resonant background, is not always 
possible. For instance, the extraction of experimental 
widths is usually model-based and relies on approxima¬ 
tions [16, 17]. As discussed in Ref. [15], many theoretical 
descriptions of collective states in the continuum [18, 19] 
employ bound-state techniques using localized wave func¬ 
tions; those have limited applicability when it comes to 
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unbound states and broad resonances in particular. 

The purpose of this work is to shed light on the notion 
of nuclear collective motion in the continuum. Ideally, 
to address this question, one should use a microscopic 
many-body approach that allows the cluster degrees of 
freedom to naturally emerge from the A-body problem 
in the open-system formulation [20-22]. The purpose of 
this work, however, is not to solve the many-body prob¬ 
lem in its full glory with the inclusion of continuum space, 
but rather to gain insight into the notion of collective res¬ 
onances by inspecting a simple case. To this end, we use 
a schematic particle-plus-core coupled-channel approach, 
based on the Berggren ensemble [23]. This model ex¬ 
plicitly contains bound states, resonances and complex- 
energy scattering states, i.e., all ingredients needed to 
describe the coupling to the one-nucleon continuum, and 
it also allows for a simple explanation in terms of famil¬ 
iar geometrical terms. In a similar context, this approach 
has been applied successfully to the description of the res¬ 
onant spectrum of the dipolar anion [24], which can be 
considered as an extreme one-electron halo. It has been 
shown that below the dissociation threshold, the motion 
of the valence electron is strongly coupled to the collec¬ 
tive rotation of the molecule, forming rotational states. 
Above the ionization threshold, however, a rapid decou¬ 
pling of the electron’s motion from the rotational motion 
of the molecule takes place, thus leading to a disappear¬ 
ance of a collective rotational band. This observation 
begs the question whether such a demise of collective ro¬ 
tation could also be expected in nuclear halo systems. 
Here, we investigate qualitatively the existence of nu¬ 
clear rotational states in the continuum of the deformed 
one-neutron halo nucleus ^^Be. 

Model and parameters - To solve a problem of a par¬ 
ticle coupled to a collective core, we apply the non- 
adiabatic coupled-channel method of Refs. [25-28] origi¬ 
nally developed in the context of deformed proton emit¬ 
ters. We assume that the positive-parity ground-state 
band of ^^Be can be viewed as a weakly-bound/unbound 
neutron coupled to a deformed core of ^°Be [29, 30]. The 
exact form of the deformed pseudo-potential represent¬ 
ing the neutron-core interaction is not essential for the 
purpose of our qualitative discussion, as long as the one- 
neutron threshold is reproduced. Here we approximate 
it by a deformed Woods-Saxon (WS) potential with a 
spherical spin-orbit term [26]. The total angular mo¬ 
mentum of the system is J = j + j^, where j = £ + s is 
the anglar momentum of the valence nucleon and is 
that of the rotor. 

In the coupled-channel (CC) formalism, eigenstates of 
the decaying nucleus ) are expanded in the basis 
of channel wave functions labeled by channel quantum 
numbers c= {£jjr)- Each channel state is given by the 
cluster radial wave function Uc{r)lr representing the rel¬ 
ative radial motion of the particle and the core, and the 


orbital-spin part \j{£,s)jr]JMj). 
The CC equations for Uc{r) are: 


2/x dr'^ 2fj,r^ 
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where is the diagonal part of the neutron-core poten¬ 
tial, Q{J^,jr) = E'^ - E^/ is the energy of the particle 
decaying from the state E'^ of the parent nucleus to the 
state E^^ of the daughter system, and V^^,{r) are the 
off-diagonal coupling terms. 

The CC equations (1) can be solved using the Berggren 
expansion method (BEM) [24, 31]. This approach has 
been benchmarked in Ref. [31] by using the traditional 
technique of direct integration of CC equations (DIM), 
which is fairly accurate for well localized states, but less 
accurate for unbound Gamow states or when a large num¬ 
ber of channels is considered. Indeed, the DIM and BEM 
results are very close for the bound ground state and low- 
lying narrow resonances, but for broader resonant states 
the BEM predicts smaller energies and widths; hence, it 
is favored by the generalized variational principle. More¬ 
over, the DIM requires a starting point for the integra¬ 
tion that is close to the final eigenenergy, otherwise it can 
sometimes exhibit a divergent behavior, and the search 
for a good starting point has to be carried out for all 
individual eigenstates of interest, while in the BEM the 
whole spectrum is obtained in one diagonalization. Fi¬ 
nally, BEM calculations are much faster, especially for 
a large number of channels. For instance, for the 7/2 
yrast state of ^^Be with ^max = 6 (22 channels), the time 
required by the DIM to solve the CC equations is over 
three orders of magnitude longer than in the BEM. 

In the present case, a Berggren basis is built for each 
partial wave {£, j), from the discrete solutions of the diag¬ 
onal part of Eq. (1) with outgoing boundary conditions. 
The off-diagonal matrix elements are calculated us¬ 
ing the exterior complex-scaling technique [32, 33]. 

The rotational structure of nuclear states can be inter¬ 
preted in terms of the internal density of the valence par¬ 
ticle in the core reference frame. Here we assume that the 
core is associated with the rigid rotor axially deformed 
around the z-axis. When expressed in the deformed ref¬ 
erence frame, the eigenstates ) can be expanded in 
the basis ), where K is the projection of the total 
angular momentum J on the symmetry axis in the core 
frame. The density pjK{r,0) is calculated as an average 
over all orientations of the density operator [24]. The 
total density is then obtained by summing up of all K- 
components. If only one AT-component is nonzero, K be¬ 
comes a good quantum number and the strong coupling 
(adiabatic) limit is strictly obeyed [25, 26, 28, 34]; the 
resulting density is referred to as the intrinsic density. 
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Coupled-channel equations were solved up to a maxi¬ 
mal radius of i?max = 30 fm and the rotation radius for the 
exterior complex-scaling i?rot = 15 fm. Complex-energy 
scattering states entering in the Berggren basis are se¬ 
lected along a contour in the complex momentum plane 
defined by three segments connecting the points ( 0 , 0 ), 
(0.2,-0.2), (0.5,0), and (2,0); allinfm”^. The contour has 
been discretized with 60 points using a Gauss-Legendre 
quadrature. Since the rotational bands in question have 
positive parity, we took partial waves with .^ = 0 , 2 ,4 ,6 in 
the Berggren ensemble, and the maximum angular mo¬ 
mentum of the core that was large enough to guarantee 
that the number of included states in the ground-state 
band of the daughter nucleus would not impact the cal¬ 
culated widths [24, 26]. In order to partly satisfy the 
Pauli principle between core and valence particles, the 
core neutron Osi /2 shell has been excluded in the con¬ 
struction of the CC basis [35]. We checked that with this 
choice of parameters, the solutions of CC equations are 
perfectly stable. 

For the energies of the core nucleus , we took the 
known experimental = 0 ^, 2 ^, and 4^ members of the 
ground-state band of ^°Be [36]. The higher-lying band 
members are approximated by means of the rigid rotor 
expression with the moment of inertia corresponding to 
the 4^ level. In the following discussion, we ignore the 
particle width of the core states with jr > 4. (Experi¬ 
mentally the 2 { state in ^°Be is particle-bound, and the 
4i level has a fairly small width of 121 keV.) 

The parameters of the deformed pseudo-potential have 
been fitted to the 1 / 2 ^ and 5/2^ members of the yrast 
band of ^^Be [36]. We have checked that these states 
collapse to the same band-head energy in the adiabatic 
limit (/ ->• oo), i.e., they are members of the same ro¬ 
tational band. It is worth noting that the location of 
higher-lying band members is very uncertain, although 
candidates for the 3/2];, 7/2];, and 9/2]; states have been 
suggested [37, 38]. The optimized parameters of the 
WS pseudo-potential are: radius Rq = 2.08 fm; diffuse¬ 
ness a = 1.1 fm; strength Vq = -59.36 MeV; spin-orbit 
strength 14o = 15.09 MeV; and quadrupole deformation 
(32 = 0.53. We have checked that the general conclusions 
of our study are independent on the precise values of 
these parameters, as long as the one-neutron threshold is 
reproduced. For instance, similar results were obtained 
using a = 0.77 fm, Rq = 2.55 fm, and (32 = 0.52. 

Results - The calculated lowest states of ^^Be are 
shown in Fig. 1. The large splitting between the favorite- 
signature band (r = exp(-i 7 r J) = -i) and unfavored band 
(r = +i) is consistent with the results of a microscopic 
multicluster model [30, 39] and large-scale shell model 
[19, 40]. It is interesting to note that Q{J,jr) < 0 for 
IJ - jr\ = 1/2 for the yrast states in ^^Be. For instance, 
the 3/2/ and 5/2/ levels of ^^Be are predicted to lie be¬ 
low the yrast 2^ state of ^°Be. This means that the i = 0 
neutron emission channel is blocked for both r = -i and 
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FIG. 1. (Color online) Calculated yrast bands of ^^Be with 
J < 13/2 with signature r = -i (squares) and r = +i (triangles) 
compared to the experimental ground-state band of ^°Be core 
(horizontal dashed lines). Some excited (yrare) states of ^^Be 
are marked by stars. 


r = +i bands. 

The rotational structure of the ground-state band of 
^^Be is revealed by looking at the weights of individual 
iF-components of valence neutron density. It turns out 
that the K = 1/2 component is dominant in most cases; 
in particular for the favored band. For the 7/2^, 11/2’*', 
and 15/2^ states, the K = 5/2 and 3/2 components dom¬ 
inate. In most cases, an appreciable degree of E-mixing 
is predicted. This suggests that a “E = 1/2” label often 
attached to this band should be taken with a grain of 
salt. 

To get further insight into the structure of the yrast 
line of ^^Be, the density pj_R'(r) of the valence neutron 
in the core frame is plotted in Fig. 2 for selected states 
and E-components. In the adiabatic limit, we checked 
that all band members have indeed the same E = 1/2 
intrinsic density shown in Fig. 2(a), whereas all E > 1/2 
components vanish. For the favored band, for which the 
E = 1/2 component dominates, the densities Pjk=i/2 for 
J = 9/2,13/2, and 17/2 are similar to that for the 5/2^ 
state shown in Fig. 2(d). This is not the case for the 
7/2^, 11/2’*', and 15/2'*' band members, which show an 
appreciable E-mixing. Indeed, as illustrated in Figs. 2(e) 
and (f), the dominant neutron distributions for the 7/2^ 
and 15/2^ states are very different. Figure 2 also shows 
that the calculated valence neutron states in ^^Be are all 
fairly well localized within the range of the neutron-core 
potential. We checked that this holds for the yrare states 
as well. 

Figure 3(a) shows the (real part) of the norms ngj = 
from various partial waves (ij) to the yrast 
band of ^^Be. It is seen that the alignment pattern of 
the valence neutron is governed by a transition from the 
Sij 2 wave, which dominates at low spins, to 1 ^ 5 / 2 , which 
governs the rotation at higher angular momenta. In the 
high-spin region, J > 7/2, the yrast line of ^^Be can be 
associated with the weak coupling of neutron’s angular 
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FIG. 2. (Color online) Partial densities pjK{r) (in 10^^ fm~^) 
with If = 1/2 for the three first ground-state band members of 
^^Be in panels (b), (c) and (d), and for selected states of the 
unfavored r = i band (in 10^^ fm~®) in panels (e) J = 7/2, if = 
5/2 and (f) J = 1512, K = 3/2. The adiabatic limit (1 ^ oo) is 
shown in panel (a). 


momentum, j = 5/2, to the angular momentum jr of the 
core resulting in the full rotational alignment of j with 
j^. Indeed, as seen in Fig. 1, the computed energies E'^ 
of the r = -i band appear close to the energies of ^°Be 
with jr = J - 5/2. The higher partial waves, such as gg/ 2 , 
appear for J > 15/2. The structure of yrare states shown 
in Fig. 3(c) is dominated by (^ 5/2 at low spins, and by 
Si /2 at higher angular momenta. 

To estimate one-neutron decay rates, we compute de¬ 
cay widths. An accurate way to compute the total 
decay width is to use the so-called current expression 
[26, 27, 41], that gives F = EcrcC?") in terms of partial 
widths from the channel wave functions. We checked that 
in every case considered the value of F obtained from the 
current expression agrees with the eigenvalue estimate 
-21ui{E^ ) and that the values of rc(r) are stable at 

^ “ 2^niax- 

The calculated one-neutron widths corresponding to 
different partial waves are shown in panels (b) and (d) of 
Fig. 3. As discussed earlier, the s-wave neutron decay is 
blocked in the ground band of ^^Be. Consequently, the 
neutron widths of yrast states are primarily governed by 
£ = 2 waves. Moreover, because of the weak coupling of 
the valence neutron to the core states, the Q-values for 
the neutron decay of the favorite band is small in the 
d^l 2 channel. As a result, the states in the favorite yrast 
band of ^^Be are predicted to have small neutron widths 
of the order of 200 keV. In the case of the unfavored 
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FIG. 3. (Color online) Contribution from different partial 
waves (£j) to the norms (top) and widths F^^ = 

(bottom) of different states of the ground-state band of ^^Be 
(left) and yrare states (right). The total widths are marked 
by stars. 


band, the neutron widths are larger, around F = 0.7 MeV. 
This result demonstrates that angular momentum align¬ 
ment can stabilize collective behavior in highly excited 
yrast states of a neutron drip-line system. The under¬ 
lying mechanism is similar to that discussed in Ref. [42] 
in the context of the cranking description of rotational 
properties of neutron drip-line nuclei. Due to the Corio¬ 
lis force, high-£ orbits responsible for angular momentum 
alignment become occupied at high spins at the expense 
of low-£ states. The latter govern halo properties and 
particle decay; the former are well localized within the 
nuclear volume because of their large centrifugar barrier. 

The non-resonant continuum does impact the pre¬ 
dicted widths of broader resonances. This is illustrated 
in Fig. 4, which displays the partial wave contributions to 
the total width for the 7/2/ state. The calculations were 
performed for £max = 6 using the flux formula in BEM, 
in both the full space and the pole space, in which the 
non-resonant contour is ignored and the basis complete¬ 
ness is broken. It is seen that removing the non-resonant 
part from the Berggren ensemble yields a width that is 
about 17% larger than the value obtained in the full ba¬ 
sis. Similar shifts can be observed for other states. 

The above discussion does not hold for excited (yrare) 
states of ^^Be with J = 1/2 and 3/2, for which the £ = 0 
channel is not blocked. As seen in Fig. 3(d), the total 
widths of those states is predicted to be around 0.4 MeV, 
and it is even totally dominated by Si /2 for J = 1/2. For 
the J = 5/2 yrare state, the £ = 0 channel is blocked again; 
hence, its width is small. 

Conclusion - We have investigated the existence of ro¬ 
tational bands embedded in the particle continuum us¬ 
ing a non-adiabatic particle-plus-core model. As case in 
point, we considered the yrast band of the one-neutron 
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FIG. 4. (Color online) Partial wave contribution to the width 
of the 7/2^ yrast state of '^^Be. Continuous and dashed lines 
correspond to the full and pole space calculations, respec¬ 
tively. The total width is indicated with dotted line in both 
cases. 


halo nucleus built upon the = 1/2^ ground state. 
We show that strong Coriolis effects result in a rotational 
alignment of the valence neutron. Owing to their large 
^-content, and the fact that low-^ channels are blocked, 
the levels forming the yrast line of ^^Be are expected to 
have small one-neutron widths: the neutron decay is pre¬ 
dicted to proceed primarily via the unfavored sequence. 
While our conclusions have been obtained by means of 
a schematic model, they nicely complement large-scale 
no-core shell model calculations [19, 40], in which a de¬ 
coupled K'^ = 1/2^ band appears in ^^Be viewed as an 11- 
body closed quantum system. Our results suggest that 
the coupling to the one-neutron continuum will play a 
minor role in the microscopic description of the yrast se¬ 
quence. This is not the case for the yrare states of ^^Be, 
for which the s-wave decay is allowed. Here, the neu¬ 
tron widths can be large and continuum coupling effects 
- appreciable. Unlike in the case of dipolar anions, where 
the valence electron density in unbound states can be ex¬ 
tremely diffused due to the shallowness of the electron- 
molecule pseudopotential, the calculated valence neutron 
states in ^^Be are fairly well localized within the range 
of the neutron-core potential. 

In summary, we demonstrated that collective states 
can exist at high excitation energy in weakly bound neu¬ 
tron drip-line nuclei such as ^^Be. The calculated decay 
widths of those states can be narrow, and the resulting 
half-lives can be long compared to the single-particle time 
scale. This justifies the use of a geometric picture in such 
cases. However, whether the notion of collective nuclear 


states applies to very broad resonances with T 112 < Tg.p. 
is still under discussion. This matter will be a subject of 
forthcoming investigations. 
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